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Percolation and Deconfinement in SU(2) Gauge Theory * 
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We show that deconfinement in SU(2) gauge theory can be described by the percolation of site-bond clusters 
of like-sign Polyakov loops. In particular, we find that in 2+1 dimensions the percolation variables coincide with 
those of the 2-dimensional Ising model. 



1. INTRODUCTION 

The critical behaviour of the Ising model can 
today be reformulated in terms of percolation the- 
ory: magnetization sets in when suitably defined 
clusters of parallel spins reach the dimensions of 
the system [0. In particular, the critical expo- 
nents for percolation then become equal to the 
Ising exponents. 

We extend this description of critical behaviour 
in terms of percolation to the deconfinement tran- 
sition in SU(2) gauge theory for two space dimen- 
sions. We show that the percolation of Polyakov 
loop clusters (taken to be suitably defined areas 
of Polyakov loops L of the same sign) leads to 
the correct deconfinement temperature and to the 
correct critical exponents for the deconfinement. 

In contrast to the conventional studies of L{T), 
the use of percolation strenght P{T) as deconfine- 
ment order parameter remains possible also in the 
presence of dynamical quarks and could thus con- 
stitute a genuine deconfinement order parameter 
for full QCD. 

2. PERCOLATION AND CRITICAL 
PHENOMENA 

The percolation problem Q is easy to formu- 
late: just place randomly pawns on a chessboard. 
Regions of adjacent pawns form clusters. Per- 
colation theory deals with the properties of these 
clusters when the chessboard is infinitely large. If 



one of the clusters spans the chessboard from one 
side to the opposite one, we say that the cluster 
percolates. 

Quantitatively, one counts how many pawns 
belong to each cluster and calculates two quanti- 
ties: 

• The average cluster size S, defined as: 



(1) 



Here rig is the number of clusters of size s 
and the sums exclude the percolating cluster; this 
number indicates how big on average the clusters 
are which do not percolate. 
• The percolation strength P, defined as: 



P = 



size of the percolating cluster 
no. of lattice sites 



(2) 
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By varying the density of our pawns, a kind of 
phase transition occurs. We pass from a phase 
of non-percolation to a phase in which one of the 
clusters percolates. The percolation strength P is 
the order parameter of this transition: it is zero 
in the non-percolation phase and is different from 
zero in the percolation phase. It is particularly in- 
teresting to study what happens near the concen- 
tration pc where for the first time a percolating 
cluster is formed. 

It turns out that P and S as functions of the 
density behave respectively like the magnetisa- 
tion M of the Ising model and its magnetic sus- 
ceptibility X as functions of the temperature T 
(Table 1). 
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Table 1 




Behaviour around the critical point 


Percolation 


Ising Model 




M(x{Tc - Tf 


for p> Pc 


for T<T, 


Soc\p-pc\~'^ 


X^\T,-T\-' 



3. DROPLET MODEL 

The analogy of percolation with second order 
thermal phase transitions leads to a natural ques- 
tion: 

Is it possible to describe these transitions 
in terms of percolation ? 

The answer to that is given by the droplet 
models which reproduces precisely the results of 
the Ising model. 

The droplet model establishes the following 
correspondence: 



spins up (down) 
spins down (up) 
spontaneous 
magnetization M 

susceptibility x 



occupied sites 
empty sites 
strength of the 
perc. cluster P 

av. cluster size S 



The clusters are basically magnetic domains 
(either with positive or with negative magneti- 
zation). When there is an infinite cluster (perco- 
lation), magnetization becomes a global property 
of the system. 

One has then to find at which temperature per- 
colation occurs and determine the corresponding 
critical exponents. But this requires one further 
conceptual feature. If the droplets are defined as 
clusters of nearest-neighbour spins of the same 
type (pure site percolation), then they do not 
lead to the correct Ising results, neither in two di- 
mensions, where the critical points coincide but 
not the exponents |3|, nor in three dimensions, 
where even the thresholds are different In 
other words, correlations in pure site percolation 
differ from thermal Ising correlations. 

To solve the problem, Coniglio and Klein 
introduced a bond probability CK = 1 — e~ 
(J is the spin-spin coupling of the Ising model, 
/3=l/kT, being k the Boltzmann constant). The 



term 2J is just the difference IS.E of energy 
that we have when we pass from a pair of next- 
neighbour spins to a one in which one of them is 
flipped. To define the new "droplets" one has to 
connect the parallel nearest-neighbour spins in a 
cluster with a bond probability CK. Since CK is 
less than one, not all the parallel spins in a cluster 
become part of the same droplet. 

With the Coniglio-Klein definition, the droplet 
model reproduces the results of the Ising model, 
both in two and in three dimensions. 



4. EXTENDING PERCOLATION 
SU(2) GAUGE THEORY 



TO 



As it is well known, considerations of symmetry 
led to the conjecture that SU(2) Gauge Theory 
and the Ising model belong to the same univer- 
sality class, that is they have the same critical 
exponents (H). This close relation inspired our 
work. 

If we take a typical SU(2) configuration at a 
certain temperature, there will be areas where L 
takes negative values, and areas where L takes 
positive values. Both the positive and the nega- 
tive 'islands' can be seen as local regions of decon- 
finement. But as long as there are finite islands 
of both signs, deconfinement remains a local phe- 
nomenon and the whole system is in the confined 
phase. When one of this islands percolates, that 
is it becomes infinite, then we can talk of decon- 
finement as a global phase of the system. 

The main point of our work is to look for a suit- 
able definition of droplets for SU(2). We expect 
our definition to be similar to the Coniglio-Klein 
one. The problem is then to find the right bond 
probability. In contrast to the Ising model, in 
which the spins can only take two values, we have 
now for L a continuous range of values. There- 
fore, if we try to build a bond probability like 
CK = 1 — e~^^^, AE is not the same at any lat- 
tice site. Because of that we used a local factor 



1 



(3) 



Here Peff is an effective coupling which was 
shown in certain cases to approximate SU(2) pure 
gauge theory as a system of nearest-neighbour in- 
teracting spins m. 
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5. RESULTS 

We performed simulations of SU(2) gauge the- 
ory in 2+1 dimensions (2 space dimensions and 1 
time dimension) using four different lattice sizes, 
642, gg2^ 128'^, 160^. We focused on the case of 
Nt — 2 lattice spacings in the time direction. 
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Figure 2. Best fit values for the exponents' ratios 
f]/]^ and calculated around the critical point. 
The vertical line indicates the /3 value where we 
got the smallest for the fit. At this point both 
P/v and take the Ising values. 



Figure 1. Average Cluster Size S of our SU(2) 
droplets for the four lattice sizes we used. The 
peaks approach the physical threshold (dotted 
line) for larger lattices. 

Figure |^ shows the average cluster size S for 
our four lattices at different (3 values. The curves 
clearly peak close to the physical transition, and 
they shift slightly to the right the bigger the lat- 
tice size is. The second step was to perform other 
simulations in the very narrow range of /3 values 
where our transition seems to occur. We per- 
formed a finite size scaling analysis of our data 
in this smaller range. Figure ^ shows the re- 
sults of our analysis. The two curves represent 
the best fit values for the ratios (3/i^ and j/i/ cal- 
culated for different beta values in the range be- 
tween f3 = 3.410 and /3 = 3.457. We found the 
best at /3 = 3.4407, where P/v and j/v take 
the Ising values 1/8 and 7/4, respectively. 

6. CONCLUSIONS 

We have shown that the deconfinement tran- 
sition for SU(2) gauge theory in 2-1-1 dimensions 
can be described in terms of Polyakov Loop per- 
colation. 
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